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. Abstract. Standard big bang nucleosynthesis (SBBN) has been re- 

CN I markably successful, and it may well be the correct and sufficient account 

r\ . of what happened. However, interest in variations from the standard pic- 

D ' ture come from two sources: First, big bang nucleosynthesis can be used 

to constrain physics of the early universe. Second, there may be some dis- 
I crepancy between predictions of SBBN and observations of abundances. 

Various alternatives to SBBN include inhomogeneous nucleosynthesis, nu- 
cleosynthesis with antimatter, and nonstandard neutrino physics. 

^ ' 

o _ 

^ ! 1. Introduction 

I The success of standard big bang nucleosynthesis (SBBN) in predicting the 

observed abundances of the light elements has led to the widespread view that 
SBBN must be correct. According to this view, any remaining disagreements 
must be due to systematic errors in observations or incorrect, or too crude, 
O I chemical evolution models. While this view may well be the right one, we should 

^ ' not be blind to other possibilities. 

. I will stay within the context of the Hot Big Bang (for an alternative, see 

Burbidge & Hoyle 1998), and discuss some models of nonstandard big bang 
nucleosynthesis (NSBBN). NSBBN scenarios range from small modifications to 
^ ' SBBN to a complete change in the decisive physical phenomena, like in the 

late-decaying massive particle scheme of Dimopoulos et al. (1988). 

Motivations for studying NSBBN go in two directions. First, the remark- 
able success of SBBN allows one to severely constrain the physics of the early 
universe. If one tries to change the conditions from the standard assumptions 
the resulting abundances of the light elements differ from the observed ones. For 
many modifications, BBN provides the strongest constraints. BBN gives also 
the strongest constraint on the single parameter of SBBN, the baryon density, 
usually given as the baryon-to-photon ratio, 

ry=^, 7710 = 101%. (1) 

Second, one may try to improve on SBBN. From time to time it has seemed that 
there might be some discrepancy between observations and SBBN, which could 
then be explained by NSBBN. In particular, there has been tension between 
D/H and Yp (see, e.g., Hata et al. 1995). To relieve this tension, either a lower 
D or a lower '^He yield has been looked for. Also one may want to relax the 
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SBBN bounds to rj. Other astronomical considerations have given motivation 
for trying to raise the upper hmit to ij. If one beUeves that the energy density 
of the universe is dominated by vacuum energy (the cosmological constant) and 
accepts the newer observations on D /H and Yp favoring somewhat larger rj within 
SBBN, this motivation largely disappears. 

There is a very large body of work on NSBBN. Extensive reviews are given 
by Malaney & Mathews (1993) and Sarkar (1996), which contain, respectively, 
over 500 and over 700 references. Here I will be able to mention only a random 
few. 

Most of the work on NSBBN can be divided into four broad classes: 

1. Inhomogeneous BBN. Usually this means inhomogeneity in the baryon-to- 
photon ratio, r/, but there are also other possibilities, like inhomogeneity 
in the neutrino chemical potentials. 

2. Nonstandard neutrino physics, e.g., additional ("sterile") neutrino flavors, 
neutrino degeneracy (asymmetry), massive i^,-, or neutrino oscillations. 

3. Late-decaying (r = 1-10^ s) massive particles, black holes, cosmic strings, 
etc. 

4. Time-varying fundamental constants. 

In the interest of time and space, I will discuss the first two classes only. 

2. Inhomogeneous Big Bang Nucleosynthesis 

The single parameter of SBBN is the baryon-to-photon ratio r], or the density 
of baryonic matter. In inhomogeneous big-bang nucleosynthesis (IBBN) one 
assumes that 77 is inhomogeneous. To get a significant effect on BBN this inho- 
mogeneity has to be large, Sri/rj ^ 1. Since the baryons make an insignificant 
contribution to the energy density at nucleosynthesis time, the total energy den- 
sity may still be essentially homogeneous. The inhomogeneity could be caused 
by, e.g., first-order phase transitions. The distance scale of this inhomogeneity 
is of crucial importance for IBBN. Without inflation, causal physics can only 
produce significant inhomogeneity at subhorizon scales (see Table 1). 

Table 1. The approximate temperature and horizon scale (in comov- 
ing units) for various events in the early universe. 



event 


T 


horizon 


EW phase transition 


100 GeV 


10"^ pc 


QCD phase transition 


150 MeV 


1 pc 


^He synthesis 


70 keV 


1 kpc 



Mechanisms connected with inflation can produce inhomogeneity at any 
scale. The isotropy of the cosmic microwave background (CMB) rules out sig- 
nificant inhomogeneity at ^ 10 Mpc scales, and it is difficult to construct an 
acceptable IBBN scenario which would explain inhomogeneity in observations. 
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In the usual IBBN models one considers a significantly smaller distance scale, so 
that while 77 is inhomogeneous during BBN, resulting in inhomogeneous abun- 
dances at first, everything gets mixed and becomes chemically homogeneous 
before or during galaxy formation. Thus the observable primordial abundances 
are homogeneous, while different from the SBBN predictions. 

The simplest version of IBBN is one where SBBN occurs with different r] 
in different parts of the universe, and the yields get mixed afterwards, so that 
one obtains the IBBN results by averaging SBBN results over the ij distribution, 
whose average we denote by f]. This kind of IBBN has a long history. Typically 
Yp goes up, ^Li goes up (down for small fj), and D goes up for large fj, and 
down for small 77, compared to SBBN with ij = f]. Leonard & Scherrer (1996) 
concluded that this way one can reduce the lower bound to r] from observations 
(in fact remove it, if arbitrary r] distributions are allowed), but the upper bound 
is essentially unchanged from SBBN, as ^Li and ^He are overproduced for larger 
77. The tension between D an ^He is worsened at the large end of the SBBN 
acceptable range. Thus this kind of modification to BBN appears undesirable. 

2.1. Small Scale Inhomogeneity and Neutron Diffusion 

The above applies to inhomogeneity with distance scales significantly larger than 
the neutron diffusion scale (~ 0.1 pc). If there is inhomogeneity at smaller scales, 
neutrons will diffuse out of the high density regions resulting in an inhomoge- 
neous n/p ratio. Especially if this results in n/p > 1 in some regions, the 
consequencies for BBN may be dramatic. This scenario (Applegate, Hogan, Sz 
Scherrer 1987) looked very exciting about ten years ago when it was noted that 
the QCD (quark-hadron) transition seemed likely to produce strong inhomo- 
geneity at just the right distance scale, and early IBBN calculations indicated a 
large reduction in Yp and increase in D/H allowing very large r/, even a critical 
density in baryons only. More detailed calculations showed that the effects were 
less dramatic, and the upper limit to i] given by D/H and Yp is raised at most by 
a factor of 2 or 3 as compared to SBBN, and this only if the inhomogeneity was 
at near the optimal distance scale (10~^ . . . 10~^ pc), and most of the baryon 
number was in the high density regions. The most severe problem for this kind 
of IBBN is ^Li overproduction. Some ^Li depletion (by a factor of 2 or 3) in Pop 
II stars is needed to allow for larger r] than in SBBN. Figure 1 is from a recent 
review of this scenario by Kainulainen, Kurki-Suonio, &: Sihvola (1999). 

Recent lattice QCD calculations favor a much smaller distance scale, al- 
though uncertainties are big enough so that the optimal distance scale cannot 
be ruled out. The distance scale from the electroweak (EW) phase transition 
must be so small that the effects on BBN cannot be large; in the best case they 
could be comparable to other small effects that have recently been included in 
accurate BBN codes. 

2.2. Regions of Antimatter 

A less-studied variant of IBBN is one where ij is allowed to have negative values, 
i.e., there are antimatter regions. This is possible in some baryogenesis scenarios 
(Dolgov 1996). Antimatter in cosmology has been reviewed by Steigman (1976). 
If the distance scale of antimatter regions is small, antimatter and matter will 
mix and annihilate in the early universe, and the presence of matter today 
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Figure 1. The ^He, D, and ^Li yields from small-scale inhomogeneous 

nucleosynthesis runs with a centrally condensed geometry, with density 
contrast R = 800 and high-density volume fraction = 0.125. The 
contours of (a) Yp, (b) log^o^^/H, and (c) log^Q ^Li are plotted as a 
function of the average baryon-to-photon ratio ij and the distance scale 
r of the inhomogeneity. The two horizontal dashed lines denote the 
horizon scale £h at the QCD (upper) and EW (lower) phase transitions. 
Prom Kainulainen et al. (1999). 
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Figure 2. The (a) ^He and (b) ^He yields function of the an- 
timatter/matter ratio R and the antimatter domain radius r. The 
distance scales are given both at T = 1 keV (in meters) and today (in 
parsecs). We plot contours of Yp and (the logarithm of) the number ra- 
tio ^Hc/H. The dotted lines show contours of the "median annihilation 
temperature", i.e., the temperature of the universe when 50% of the 
antimatter has annihilated. Typically the annihilation is complete at 
a temperature lower than this by about a factor of 3. The dot-dashed 
line gives the upper limit to R from CMB spectrum distortion. This 
plot is for r/io = 6. Prom Kurki-Suonio & Sihvola (1999). 

implies that there was initially more matter than antimatter. If the distance 
scale is large, so that antimatter regions will survive till present, observational 
constraints require either the amount of antimatter to be very small, or the 
distance scale to be very large, comparable to the present horizon or larger 
(Cohen, De Riijula, & Glashow 1998), so that the case of large regions is not of 
interest for BBN. 

The smaller the antimatter regions are, the earlier they annihilate. Rehm 
& Jedamzik (1998) considered annihilation immediately before nucleosynthesis. 
Kurki-Suonio & Sihvola (1999) extended these results to larger distance scales 
where annihilation occurs during or after nucleosynthesis (see Figure 2). So far 
the focus has been on obtaining upper limits to the amount of antimatter at 
various scales in the early universe, but clearly there is also potential for ob- 
taining acceptable abundances with nonstandard values of r], although probably 
only with fine-tuned model parameters. 



3. Neutrinos and Big Bang Nucleosynthesis 

Neutrinos affect BBN in two ways, through the energy density effect and the Vg 
effect. The most significant effect is on Yp in both cases. 

The energy density in neutrinos affects the expansion rate of the universe. 
The simplest way to increase the energy density of the early universe from the 
standard model is to have additional particle species (sterile neutrinos or other 
hypothetical particles). The custom is to parametrize this by an "effective num- 
ber of neutrino species". The standard case is Nj^ = 3. We now know that 
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there are only three "active" neutrino species, so any additional species must be 
"sterile" neutrinos or other very weakly interacting particles. A higher energy 
density means faster expansion. This leads to n/p freezeout at a higher temper- 
ature, leaving more neutrons, and resulting in a higher '^Hc yield. The D yield 
is also increased, so an increased energy density is disfavored by BBN, and one 
gets an upper hmit, e.g., iVj, < 3.2 (Buries et al. 1999) or A^^j, ^ 4 (Lisi, Sarkar, 
& Villantc 1999), depending on what observational constraints one uses. 

Electron neutrinos affect the weak n p reactions directly. More i^e leads 
to fewer neutrons and thus to less ^He (and everything else), whereas more i?e 
leads to more neutrons and more ^He. 

3.1. Neutrino Degeneracy 

In SBBN one assumes that the neutrino asymmetry (difference between the 
number of neutrinos and antineutrinos). 



which is related to the neutrino chemical potential /i^, or the degeneracy param- 
eter = fJ-u/T, is small, <C 1. This seems natural, since the comparable baryon 
asymmetry i] is small. However, the neutrino background is unobservable, so 
we cannot rule out a large neutrino asymmetry. A larger asymmetry always 
means a larger neutrino energy density, raising N^. To have a significant effect 
on BBN, we must have |.^|, \Li,\ ^ 0.1. There is a separate contribution from 
each neutrino flavor. Thus there are three indepedent degeneracy parameters, 
^e, ^11, and ^r- The energy density effect is the same for all three flavors, and 
depends only on |^|. The electron neutrino effect depends only on ^g, but is 
much stronger, and the direction of the effect depends on the sign. 

There arc two possible scenarios for affecting BBN. If is comparable 
in magnitude to and ^r, or larger, one can forget the other two in first 
approximation. One can then adjust to dial in the desired value of Yp. The 
other elements are hardly affected. A less natural scenario is one where the 
asymmetries in the other two neutrino flavors are much larger, and the energy 
density and z^e effects are balanced against each other to keep Yp in the acceptable 
range. This way one can have a significant effect on the other abundances 
and raise the acceptable range for rj. This second scenario is constrained by 
structure formation, since the large neutrino energy density means that the 
matter /radiation equality and thus the beginning of structure formation occurs 
later. Kang & Steigman (1992) used a generous lower limit for matter/radiation 
equality, Zcq > 10^ to widen the SBBN acceptable range from r/io = 2.8-4.7 to 
mo = 2.8-19. 

3.2. Inhomogeneous Neutrino Degeneracy 

The different results from high-z D/H measurements (Tytler, Fan, & Buries 
1996; Webb et al. 1997) raised the question whether there might be a large- 
scale inhomogeneity in primordial abundances. This is very difficult to achieve, 
since the extreme isotropy of the CMB rules out any significant large-scale inho- 
mogeneity in r] or the energy density. Dolgov &; Pagel (1999) have come up with 
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a way of getting around this constraint. In their model the asymmetries of the 
different neutrino flavors are inhomogeneous but balanced with each other so 
that they add up to a homogeneous total energy density. The inhomogeneous 
is then responsible for the inhomogeneous primordial abundances through the 
effect. They suggest that an Afflect-Dine type scenario of generation of leptonic 
charge asymmetry, respecting the symmetry between different lepton families, 
could be responsible for creating a domain structure, where the neutrino asym- 
metries would have the same three values but interchanged with respect to e, /i 
and r. To achieve a significant D/H inhomogcncity, a huge Yp inhomogcncity 
has to be allowed. But since there are no high- 2: Yp determinations, this can- 
not be used to rule out their model. Table 2 shows an example of what kind 
of abundances we could have in such a domain structure. The first line would 
correspond to our local domain; from the other domains we would have only 
D/H observations. 

Table 2. Abundances of light elements for r/10 = 5 and nonzero values 
of all three chemical potentials. One example from Dolgov &; Pagel 
(1999). 









D/H 


Yp 


'Li/R 




0.1 


-1 


1 


3.8 X 10"^ 


0.23 


2.5 X 10" 


-iU 


-1 


0.1 


1 


9.2 X lO"'"^ 


0.55 


4.5 X 10" 


10 


1 


-1 


0.1 


2.8 X 10""' 


0.08 


1.1 X 10" 


-10 



3.3. Decay of a Massive Tau Neutrino 

If the rest mass of a neutrino species is much larger than 100 MeV, then it is 
becoming nonrelativistic before nucleosynthesis and its contribution to the en- 
ergy density is different from the standard zero-mass case. The laboratory limits 
for the neutrino masses leave this as a possibility for Ur- Above the neutrino 
decoupling temperature, T ~ 3 MeV, a massive neutrino species contributes less 
energy density, because of neutrino-antineutrino annihilation, but after neutrino 
decoupling the annihilation ceases and the rest mass then contributes extra en- 
ergy density. Neutrinos this heavy must decay to avoid contributing too much to 
the present energy density. The decay time and mode are of crucial importance 
to BBN. If the decay time is very short, then the contribution to Ni, will be less 
than one. The most interesting case is the one where Vr decays into ly^ (and a 
scalar particle), since then the ly^ effect could cause a significant reduction in Yp. 

These calculations are difficult since the decisive effects occur near the neu- 
trino decoupling temperature, so thermal equilibrium is not maintained and the 
neutrino spectra arc distorted. The recent results by Hannestad (1998) and 
Dolgov et al. (1999) are in disagreement with each other. Hannestad gets the 
maximum reduction of Yp, from the SEEN result Yp = 0.239 to < 0.20, for 
ly-r mass nii, = 0.2-0.5 MeV and lifetime r < 100 s. According to E)olgov et al., 
the maximum reduction is less, to 1^ ~ 0.21, and occurs for larger masses, rui^ 
= 2-3 MeV, and requires a shorter lifetime r < 1 s. 

The most natural explanation of the SupcrKamiokandc (1998) result on 
atmospheric neutrinos is Vp, — > Vr oscillation. Then Vt cannot be heavy and 



8 



Kurki-Suonio 



its mass will not affect BBN significantly. To allow the above scenario, the 
atmospheric neutrino oscillations would have to be into a sterile neutrino species, 
i^ii ^ ^sj instead (Kainulainen et al. 1999). 

3.4. Neutrino Oscillations 

Observations of solar neutrinos and atmospheric neutrinos (SuperKamiokande 

1998) as well as the LSND (1998) accelerator experiment see different amounts of 
the different neutrino flavors than predicted by the Standard Model. This can be 
explained by neutrino oscillations. This is a quantum-mechanical phenomenon 
where the flavor (z^g) ^^i^ ^t) content of the neutrino varies periodically. This 
requires nonzero neutrino masses and the effect is determined by the difference 
in mass-squared, Am^, and the "mixing angle". 

All three (solar, atmospheric, and LSND) "neutrino problems" cannot be 
simultaneously explained by oscillations among three flavors, but require at least 
a fourth flavor, Vg, which must be "sterile", i.e., much more weakly interacting 
than the three known "active" flavors, in order not to violate the limit Ni, ~ 3 
from decay width (Particle Data Group 1998). A sterile neutrino would also 
be useful for supernova nucleosynthesis (Peltoniemi 1996; Caldwell, Fuller, k, 
Qian 1999). 

The LSND results are controversial, so the other viewpoint is to ignore them 
until they are confirmed by independent experiments, in which case the solar 
and atmospheric neutrino problems can be explained just with the three active 
neutrinos. 

Oscillations among (light, non-degenerate, i.e., ^ = 0) active neutrinos do 
not affect BBN, since they all have equal abundances. If the sterile neutrino 
exists, it would have thermally decoupled from the other neutrinos very early, 
much before BBN, so that its contribution to Ny would be ^ 1. Active-sterile 
neutrino oscillations before BBN would then lead to production of Vg , increasing 
Ni^ (Enqvist, Kainulainen, &; Thomson 1992), which from the BBN point of 
view is undesirable. The situation is more complicated, however. The oscillation 
depends on the background temperature, and at a certain temperature there is 
a resonance. This resonance temperature depends on the neutrino energy, so 
as the temperature falls, the resonance sweeps through the neutrino spectrum. 
If there is a small pre-existing asymmetry (this will be the case, since thermal 
fluctuations suffice), the rates of neutrino and antineutrino oscillation will be 
different. Resonant active sterile neutrino oscillations will then lead to a growth 
of the neutrino asymmetry by a large factor (Barbieri & Dolgov 1991; Foot & 
Volkas 1995; Shi 1996; Enqvist, Kainulainen, k Sorri 1999; Di Bari k Foot 2000). 
This may generate a large enough electron neutrino asymmetry to affect BBN 
(Bell, Foot, k Volkas 1998; Kirilova k Chizhov 1998; Shi, Fuller, k Abazajian 

1999) . 

Depending on the oscillation parameters, the asymmetry may either just 
grow or oscillate between positive and negative values, so that the final sign 
of the asymmetry becomes unpredictable. To calculate the effect on BBN is 
complicated, since the resulting distortion of the spectrum is also important 
for BBN, and the process happens near the neutrino decoupling temperature. 
There are two schemes to generate a large asymmetry, either directly via 
Ve. <-> i^s oscillations or indirectly via Vn{r) ^ and f//(r) ^ oscillations. 
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Figure 3. The impact on the primordial ^He abundance Y if an asym- 
metry in i^eZ^e is generated by a resonant z/e I's mixing (left) or by 
the indirect neutrino mixing scheme (right), as a function of Am?. The 
baryon-to-photon ratio is set to r] = 5.1 x 10"^''. From Shi et al. (1999). 



This scenario is under active study and there is much controversy among 
the different research groups. In Fig. 3 we show results obtained by Shi et al. 
(1999). The maximal effect on Yp seems to be at the ±0.01 level. 



4. Conclusions 

At present, no NSBBN scenario appears as convincing as SBBN, which is the 
simplest of all. Often the real world has turned out to be more complicated in 
the end than first assumed, but for the early universe a simple picture has been 
very successful. However, it is healthy to keep in mind the possibility that SBBN 
might not be the full story, and that any discrepancies between observations and 
SBBN might actually be telling us something important about the early universe 
or particle physics. 

Acknowledgments. I thank K. Kainulainen, A. Kalliomaki, J. Peltoniemi, 
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